Evolution of the potential-energy surface of amorphous silicon 
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The link between the energy surface of bulk systems and their dynamical properties is gener- 
ally difficult to establish. Using the activation-relaxation technique (ART nouveau), we follow the 
change in the barrier distribution of a model of amorphous silicon as a function of the degree of 
relaxation. We find that while the barrier-height distribution, calculated from the initial minimum, 
is a unique function that depends only on the level of distribution, the reverse-barrier height dis- 
tribution, calculated from the final state, is independent of the relaxation, following a different 
function. Moreover, the resulting gained or released energy distribution is a simple convolution 
of these two distributions indicating that the activation and relaxation parts of a the elementary 
relaxation mechanism are completely independent. This characterized energy landscape can be used 
to explain nano-calorimetry measurements. 

PACS numbers: 10.4 



The concept of energy landscape has been used exten- 
sively in the last decade to characterize the properties 
of complex materials. In finite systems, such as clusters 
and proteins, the classification of energy minima and en- 
ergy barriers has shown that it is possible to understand, 
using single-dimension representations such as the dis- 
connectivity graph, the fundamental origin of cluster dy- 
namics and protein folding [JJ. For bulk systems, where 
the relations are more complex, the energy landscape pic- 
ture has helped further our qualitative understanding of 
glassy dynamics and of the evolution of supercooled liq- 
uids through their inherent structures [JJ. Similarly, ap- 
plied to amorphous semiconductors, it has revealed an 
unexpected simplicity that suggests universality [2J-0] . In 
spite of these contributions, many questions remain re- 
garding the evolution of the local structure of the energy 
landscape itself as a function of relaxation, an evolution 
that determines the dynamics of disordered materials 
away from equilibrium [5J. For example, following mea- 
surements of heat released during the annealing of dam- 
age generated by low-energy ion implantation in amor- 
phous silicon, it was suggested that the activation energy 
and the amount of heat released are uncorrelated [51 [3] , 
microscopic details were still missing, however, as well as 
a general picture of the phenomenon. Results presented 
here tie experiments and simulations and provide a clear 
link between the structure of the energy landscape and 
experimental observations, supporting the importance of 
this concept. 

More precisely, we focus on amorphous silicon, a model 
system studied extensively, both experimentally and the- 
oretically, over the years (HO [7|. Using ART nouveau [5], 
a saddle-point finding method, we characterize the evo- 
lution of the local energy landscape, defined by the dis- 
tribution of transition states and adjacent energy min- 
ima around a local configuration, as a function of re- 
laxation. The resulting theoretical picture can be used 



to understand and explain recent nano-calorimctric mea- 
surements on ion-implanted a-Si samples [5J [3] . 

Our a-Si model is relaxed using the Stillinger- Weber 
potential with parameters adjusted specifically to recover 
the structural and vibrational properties of the amor- 
phous system [5]. While this potential has a number of 
shortcomings, its structural evolution as a function of 
relaxation is in good qualitative agreement with experi- 
ment as discussed below. Moreover, it is fast enough to 
allow for an extensive search of transition states, neces- 
sary to construct reliable energy distributions. Starting 
from a 4000-atom random configuration in a cubic box 
with periodic-boundary conditions, the system is relaxed 
using ART nouveau: beginning in a local energy min- 
imum, the system is brought to an adjacent transition 
state and relaxed into a new minimum. The move is ac- 
cepted using a Metropolis criterion with a temperature 
of 0.25 eV, following Ref. 010. A total of 95 445 events 
is necessary to relax the model from an initial energy per 
atom of -2.89 eV/atom going down to -3.04 eV/atom (see 
Fig. [I] in inset). Figure [I] shows the correlation between 
the energy per atom and the bond-defect concentration 
(mostly three-fold and five- fold coordinated Si). Follow- 
ing this relaxation, we find a linear relation between the 
width of the TO peak, calculated using the relation of 
Ref. [TO] , and the configurational energy, with a slope 
of 0.989 (KJ/mol)/deg, in good agreement with experi- 
ment [51 111]. 

To characterize the evolution of the local energy land- 
scape as a function of the degree of relaxation, we explore 
extensively the transition state distributions around four 
configurations with decreasing internal stress. The en- 
ergy per atom for these configurations is indicated in 
Table [TJ From each of these four generic configurations, 
we generate 100 000 events (exit pathways), providing 
a detailed picture of the local energy landscape associ- 
ated with various degrees of relaxation. In each case, we 



2 



TABLE I: Properties of the four sampled 4000-atom a-Si con- 
figurations including energy per atom and the fitting parame- 
ters to Equation [T]for the forward-barrier (FB) distribution. 
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C3 


C4 


Energy /atom 


-2.980 


-2.996 


-3.019 


-3.042 


(E FB ) (eV) 


1.919 


2.190 


2.430 


2.587 


ofb (eV) 
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FIG. 1: Fraction of coordination defects (mostly 3-fold and 
5-fold coordinated atoms) as a function of the configurational 
energy per atom (in eV) for the 4000-atom model of a-Si. The 
large circles indicate the location of the four configurations 
selected for further analysis. Inset: energy per atom as a 
function of event during the preparation of the configurations. 



remove the duplicate events leaving about 12 000 differ- 
ent activated pathways from the initial local minimum 
to a nearby one for each of these four configurations, 
i.e. about 12 000 different transition points and adja- 
cent minima. The nature of dominant diffusion and re- 
laxation mechanisms in a-Si have already be discussed 
elsewhere [2J [H [T2] and we focus here on the distribu- 
tion of forward and reverse barrier heights (measured, 
respectively, from the transition state to the initial and 
final state) and energy asymmetry (given by the energy 
difference between the final and initial state), for each 
level of relaxation. From previous analysis, we know that 
although the total number of events in a-Si is between 
30 and 60 per atom, the overall distribution converges 
relatively quickly ^ and we observe only negligible dif- 
ferences in the shape of the various energy distributions 
by including the events generated after 50 000 or 100 000 
ART moves. 

From transition state theory, we know that the kinetics 
of activation-dominated systems is controlled by the the 
rate of escape from the initial energy minimum. Figure 
[2]ja) shows the distribution of activation barriers from the 



initial minimum to the transition state (forward barrier 
or FB) for configurations CI to C4. We observe that 
the distribution is continuous and bounded in agreement 
with previous simulations [H H] but also that it shifts 
towards higher energy as the configuration are relaxed. 
As shown in Fig. |3ja), however, the four curves are well 
represented by the same function given by a Gaussian 
multiplied by the energy: 



G FB (E) = AEcxp<^ - 



{E-(E%)f 



(1) 



where, as is shown in Table |TJ the peak position con- 
trolled by E r F e B varies with the degree of relaxation but 
the overall scale A and the width of the distribution, 
given by (Tfb, are constant. Remarkably, the shift in EJfg 
is significantly larger than the drop in the average strain 
per atom as the configuration is relaxed: E r p B shifts by 
0.67 eV from CI to C4 as the energy per atom decreases 
by only 0.06 eV. This result is consistent with the grow- 
ing super-Arrhenius relaxation time scale observed in 
strong glasses and generally described by the empirically- 
derived Vogel-Fulcher law [13]. However, while this rela- 
tion postulates a constant energy-barrier and an unusual 
temperature relation, we see here that the reverse is hap- 
pening: it is the barrier height that increases with relax- 
ation. The disproportionate shift in the forward barrier 
distribution as a function of relaxation can be can be ex- 
plained by previous observations that the defected and 
strained structures play a major role in the evolution of 
amorphous systems as diffusion or rearrangement in per- 
fectly coordinated region of the amorphous network tend 
to have an energy cost similar to the of the crystal [21 [7] . 

Even though a system's kinetics is conditioned by the 
barrier height between minima, the long time evolution 
of a system is also determined by the overall event en- 
ergy asymmetry. For example, if system jumps onto a 
very energetic state, with a high energy asymmetry with 
respect to the initial minimum, it will come back imme- 
diately onto the initial state irrespective of the height of 
the forward energy barrier. In effect the presence of these 
highly asymmetric states surrounding a local minimum 
only slows down the relaxation. To characterize this phe- 
nomenon, we plot the reverse energy barrier distribution 
(Fig. [2^b)), i.e., the energy barrier height measured from 
the final minimum. As had been observed before, we note 
that the forward and reverse energy-barrier distributions 
are qualitatively different [4] : instead of a Gaussian, we 
observe an exponential distribution with a large propor- 
tion of irrelevant states and very few low-energy minima 
contributing to the evolution of a-Si into a more stable 
local state (also observed for Lennard- Jones and a-Si in 
a slightly different set-up in Ref. [3j [4] ) . Contrary to the 
behavior of the forward energy distribution, we find that 
the reverse energy barrier distribution is essentially in- 
dependent of the relaxation level, with the exception of 
region of very low energy, zero eV where the identifica- 
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tion of recurring events is sometimes tricky. Neglecting 
this region, we see that all backward barrier distributions 
are represented by the same exponential, exp(— AE/E ), 
where E = 0.60 eV (see Fig. [jjfb)). 
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FIG. 2: (a) Forward and (b) reverse energy-barrier distribu- 
tion, and (c) asymmetry energy distribution for configurations 
CI to C4. 

These two distributions, forward and backward, can- 
not be measured directly experimentally. However, the 
heat released during relaxation is directly associated with 
the energy asymmetry, i.e., the energy difference between 
the initial and final minima. Fig. [2jc) shows the en- 
ergy asymmetry distribution surrounding configurations 
CI to C4, a structure that changes with the degree of re- 
laxation. While the shape of the distribution looks com- 
plex, it is a simple convolution of the two independent 
forward and backward energy barrier distributions, as is 
shown in Fig. |3^c) . Not only do these two distribu- 
tions behave differently with the relaxation, these turn 
out to be independent: the height of a barrier does not 
provide any information regarding the energy of the final 
minimum. 

The energy landscape picture developed here can be 
used to explain recent differential nano-calorimetric mea- 
surements made on implanted o-Si [5] [5] , offering an ex- 
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FIG. 3: Collapsed (a) Forward and (b) reverse energy barrier 
distributions, following the equations given in the text, (c) 
Asymmetry distribution for configuration C2 and the corre- 
sponding distribution generated by convoluting the respective 
forward and reverse barrier distributions. 



perimental check on these results. Samples are implanted 
at low temperature and then heated, while heat released 
is measured, revealing a number of features: First, the 
heat flux released as the sample is heated is relatively 
uniform as a function of temperature, for temperatures 
ranging from 118 to 775 K; second, the stored-strain, 
measured as the heat rate, saturates as a function of im- 
plantation fluence; third, for low fiuences, the heat-rate 
follows nearly the same curve, irrespective of the implan- 
tation temperature. For fiuences near saturation, the 
heat release amplitude increases with decreasing implan- 
tation temperature. These results were explained by the 
presence of a continuous distribution of states where each 
annealing event releases an amount of heat independent 
of its activation energy. Dynamic annealing was consid- 
ered responsible for the increasing signal with decreas- 
ing temperature near saturation: some configurations are 
only accessible if underlying, low activation energy con- 
figurations survived dynamic annealing. 

The saturation of stored-strain can be readily ex- 
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plained within the context of energy landscape: We have 
seen that the barrier-height distribution shifts to lower 
energy as the strain level increases. As this shifts occurs, 
the number of barriers within k^T increases reaching a 
point, for a given temperature, at which relaxation takes 
place faster than the strain increase, leading to the satu- 
ration. 



0.4 Si7nm 2 C3 




T(K) 

FIG. 4: Heat flux as a function of temperature for Configu- 
rations CI to C4 (lines) compared with the heat flux released 
during relaxation of ion-damaged a-Si samples (experimental 
data from Karmouch et al. [5]). 



Efb) and most barriers crossed are of the order of a few 
k^T. However, since most of the weight in the reverse- 
barrier distribution is concentrated near zero, the num- 
ber of crossing leading to heat released does not change 
significantly, leading to the experimentally observed uni- 
form increase as a function of temperature. We can see in 
Fig. |4j however, that the exponential fall in reverse bar- 
rier flattens the heat rate curve for well-relaxed sampled, 
corresponding to very low fluence, in agreement with ex- 
periments. 

In conclusion, a detailed analysis of the evolution of 
the energy landscape for model a-Si shows features that 
provide a clear picture of defect relaxation in support of 
recent differential nano-calorimetric measurements. Be- 
cause the forward and reverse energy barrier distributions 
are independent, heat released after ion-implantation is 
proportional to the number of barrier crossed, which in- 
creases linearly with the temperature, for highly disor- 
dered systems, and flattens out for more relaxed con- 
figurations, where the impact of the exponential reverse 
barrier distribution becomes more important. Moreover, 
the rapid shift in the forward-energy barrier distribution 
as a function of relaxation also explains the saturation 
in stored strain at high fluence. While it remains to be 
confirmed numerically, the picture developed here should 
apply to covalent glasses in general. 



Understanding the behavior of the released heat-flux 
requires more analysis. For this, we suppose that since 
the energy-barrier distribution shifts rapidly to higher en- 
ergy as the system is relaxed, low-energy barriers present 
in the initial configuration, before heating is started, are 
not replaced: once a barrier is "used" it do not reappear 
at once because heating is relaxing the system. We do 
not have to worry about the change in shape at high- 
energy because these events have an exponentially small 
probability of occurring during the fast heating. Because 
the reverse distribution is independent of the relaxation 
state, we also pose that it remains constant during the 
heating process. 

Starting from one initial minimum, we assign a 
temperature-dependent Boltzmann probability of cross- 
ing each barrier, posing a constant attempt fre- 
quency [12]. Using a constant heating rate, we select 
forward-energy barriers with the Bolztmann-weighted 
probability times the probability distribution and the re- 
verse barrier with the exponential distribution. An event 
is accepted only if the final state has an energy lower than 
the initial state, since high-energy states will be unstable 
and will jump back immediately to the original minimum. 
The resulting heat flux, shown in Fig. [4] compares very 
well with the experimental results allowing us to relate 
this behavior to two factors: first, the number of barri- 
ers being crossed as a function of temperature increases 
linearly only (because T is very small with respect to 
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